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Abstract 

The two- and all-terminal reliabilities of the Brecht-Colbourn ladder and the gen- 
eralized fan have been calculated exactly for arbitrary size as well as arbitrary 
individual edge and node reliabilities, using transfer matrices of dimension four at 
most. While the all-terminal reliabilities of these graphs are identical, the special 
case of identical edge (p) and node (p) reliabilities shows that their two-terminal 
reliabilities are quite distinct, as demonstrated by their generating functions and 

the locations of the zeros of the reliability polynomials, which undergo structural 
1 

transitions at p = —. 
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1 Introduction 

Network reliability has long been a practical issue since the pioneering work of 
Moore and Shannon [25], and will remain addressed by reliability engineers, 
statistical physicists and applied mathematicians for years, since networks are 
pervading our everyday life. Not surprisingly, the study of network reliability 
has led to a huge body of literature, which includes excellent textbooks and 
surveys [61171 26 .59,60]. In the following, we consider a probabilistic approach 
of reliability, in which the network is represented by an undirected graph 
G = (V,E), where V is a set of nodes (also called vertices), and E is a set 
of undirected edges (or links). Each element of V and E has a probability 
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p n or p e to operate correctly; failures of constituents are assumed to occur 
at random, and to be statistically independent events (this restriction may 
be relaxed, eventually). Among the different measures of reliability, one often 
considers the /c-terminal reliability, i.e., the probability that a given subset K 
of k nodes (K C E) are connected. The most common instances are the all- 
terminal reliability Rel^ (K = E) and the two-terminal reliability Rel2(s — > t), 
which deals with a particular connection between a source s and a terminal 
destination t. Both of them are affine functions of each p n and p e . 

The sheer number of possible system states, namely 2' s ' + ' v '', clearly precludes 
the use of an "enumeration of states" strategy for realistic networks, and shows 
that the final expression may be extremely cumbersome. Consequently, most 
studies have assumed graphs with perfect nodes (p n = 1) and edges of identical 
reliability p; all reliabilities are then expressed as a polynomial in p, called the 
reliability polynomial. Radio broadcast networks [35] have also been described 
by networks with perfectly reliable edges but imperfect nodes; in this context, 
one speaks of residual connectedness [12] . It was shown early on [26] that the 
calculation of /c-terminal reliability is T^P-hard, even after the simplifying and 
restricting assumptions that (i) the graph is planar (ii) all nodes are perfectly 
reliable (iii) all edges have the same reliability p. 

The difficulty of the problem has stimulated many approaches : partitioning 
techniques [3"U1T77] . sum of disjoint products [2T51I38|I53|I63"] . graph simplifica- 
tions (series-parallel reductions [16], triangle-star (also called delta- wye) trans- 
formations [20113 ll!34f 551172] . factoring [12]), determination of various lower and 
upper bounds to reliability polynomials [6|13f 14lllTf2 2 26 2 7|52f58j . Monte- 
Carlo simulations [3"3|4Uf47j . genetic [23J and ordered binary decision diagram 
(OBDD) algorithms [4"3l54f75lTT6] . The reliability polynomial has also been 
studied [51] with the aim of deriving some useful and hopefully general infor- 
mation from the structure of its coefficients [2TT28] . or the location of its zeros 
in the complex plane [T6j . 

The Tutte polynomial T(G, x, y) of a graph G has also been shown to be equiv- 
alent to the Potts model partition function of the g-state Potts model [5"TfT3"] . 
Calculations for various recursive families of graphs G quickly followed |19j . 
The all-terminal reliability polynomial Ke\A{G,p) of graphs is deduced from 
T(G, 1, jzr)- Royle an d Sokal [56] proved that the Brown- Colbourn conjecture 
[T6] on the location of the zeros of Ke\A_{G,p), while valid for series-parallel re- 
ducible graphs, does not hold for a few families of graphs. While these results 
are extremely valuable to better understand a few properties of graphs and 
all-terminal reliability polynomials, they still assume that nodes are perfect 
and that edges have the same reliability. 

In recent years, the growth of Internet traffic has called for a better evaluation 
of the reliability of connections in — among others, optical — networks. This, 
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of course, strongly depends on the connection under consideration. Actual 
failure rates and maintenance data show that a proper evaluation of two- 
terminal reliabilities must put node and edge equipments on an equal footing, 
i.e., both edge (fiber links, optical amplifiers) and node (optical cross-connects, 
routers) failures must be taken into account. The possibility of node failure 
has been considered in early papers [T1l32l36j . to quote but a few. Adapta- 
tion of algorithms to include imperfect nodes has been addressed, sometimes 
controversially |48. 4~TIT0l71|76j ; the two- variable approach for bounds to the 
reliability polynomial, by Bulka and Dugan [TH] and Chen and He [22], is also 
worth mentioning. In order to be realistic, different edge reliabilities should 
be used: for instance, the failure rate of optical fiber links is likely to increase 
with their length. 




Fig. 1. The Brecht-Colbourn ladder. 




T 

Fig. 2. The generalized fan. 



In this work, we give the exact two-terminal reliability of two well-known 
graphs, namely the Brecht-Colbourn ladder [T3lll4f52j and the generalized fan 
[3l4U| . which have the same all-terminal reliability. Both networks have been 
studied in the literature, especially the first one, as a case study for different 
lower bounds of the two-terminal reliability [2] , but with a limited number of 
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nodes. Here, the number of nodes will be arbitrary. Our solution fully exploits 
the recursive nature of the graphs and the triangle-star transformation with 
unreliable nodes |34j . The final expressions are products of (at most) 4x4 
transfer matrices, in which the arbitrary reliability of each node and edge 
appears explicitly. In the case of identical edge (p) and node (p) reliabilities, 
we provide the generating functions of the two-terminal reliabilities, leading 
to a very simple expression for the generalized fan. Another byproduct is the 
asymptotic locations of the zeros of the two-terminal reliability polynomials 
for the Brecht-Colbourn ladder and the generalized fan, which exhibit very 

different structures and undergo a structural transition at p = -. 



Our aim is (i) to give a detailed derivation of the final results, so that people 
involved in reliability studies can readily use an easy-to-implement formula (ii) 
compare the exact solution with previous results (iii) exhibit the structural 
changes undergone by the location of the zeros of reliability polynomials as a 
function of node reliability (iv) emphasize anew the importance of algebraic 
structures of the underlying graphs in the determination of their associated 
polynomials [T0f59] . 




Our paper is organized as follows : In Section [21 we briefly recall the formu- 
lae for the triangle-star transformation for unreliable nodes. In Section [21 we 
define the notations for the different edge and node reliabilities and detail the 
derivation of the main results (eqs. (fT2l) and f[T4l for the Brecht-Colbourn lad- 
der, eqs. (}23l) and f[25|) for the generalized fan). Section 0] is devoted to the case 
where all edges and nodes have identical reliabilities p and p, respectively, the 
size of the network appearing simply as an integer n. We give the analytical 
solution of the two-terminal reliability RekjQo, p;n) for the Brecht-Colbourn 
ladder and the generalized fan, and the associated generating functions, which 
encode all the necessary information in a beautifully simple, compact form. 
The asymptotic power-law or constant behaviors are given when n — ► oo. 
Prompted by the nearly universal character of the Brown-Colbourn conjec- 
ture [Ej, we address in Section \5\ the location of zeros of the two-terminal 
reliability polynomials, and show that their structures are very different even 
though their internal structure is similar. For the sake of completeness, we 
derive in Section [6] the common all-terminal reliability for both networks, for 
arbitrary values of edge reliabilities. Finally, we conclude by indicating several 
directions in which the present results may be further extended, so that, for 
instance, a catalog of exactly solvable networks — in terms of reliability - 
may be given rapidly [69]. Such a catalog of elementary bricks could be useful 
for a new and improved set of bounds or benchmarks for alternative methods 
in the general case. 
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2 Triangle-star transformation for unreliable nodes 



The triangle-star — also called delta- wye or A — Y — transformation has been 
used many times to simplify calculations of network reliability |20f26ll31f34f42ll55f72] 
It has mostly been applied in a perfect nodes context, to provide upper and 
lower bounds to the exact reliability. Here, we exploit this transformation to 
the full in the case of imperfect nodes in order to obtain exact results. Since it 
plays a crucial part of the derivation, we give the formulae derived by Gadani 




Fig. 3. Triangle-star transformation for unreliable nodes. A, B, and C are the node 
reliabilities, with a, b, and c the edge reliabilities of the initial network, and pa, Pb, 
pc, and O, those of the transformed network. 



Let us consider three particular nodes A, B, and C of a network, connected 
as shown in the left part of Fig. [3] (the reliability of the nodes are given by 
the same uppercase variables, so as to prevent an unnecessary multiplication 
of notations). The reliability of the edge connecting A and B is given by (low- 
ercase) c, with similar notation for the remaining edges of the triangle. The 
aim of the triangle-star transformation is to replace the triangle by a star, 
which is possible by the addition of a new unreliable node O and three new 
edges connecting O to A, B, and C, with reliabilities pa, Pb, and Pc, respec- 
tively. Both networks are equivalent provided that the following compatibility 
relations hold [M] 



PaO pc = b + acB — abc B, (1) 

PaOpb = c + abC — abcC, (2) 

Pb Opc = a + bcA — a bcA, (3) 

Pa OpsPc = a b + be + a c — 2 a be. (4) 

Note that the first three equalities correspond to the probability that the 
two nodes under consideration are connected, while the last one gives the 
probability that the three nodes are connected. A word of caution — already 
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given by Gadani — is worth mentioning in the case of successive triangle-star 
transformation: the triangles should have no common edge or node. 



3 Derivation of the main results 

3.1 Brecht-Colbourn ladder 

Let us first name the different edge and node reliabilities of the Brecht- 
Colbourn ladder, and detail how we can use the triangle-star transformation 
of the preceding section. Assuming that S is always the source node, and 
using node reliabilities Si, we note a n the reliability of the edge (S n -2,S n ), 
and b n that of (S n -i, S n ), as shown at the top of Fig. HI 



The application of the triangle-star transformation to Fig. H] leads to 



Pi Po P2 = b n + a n b n -i S n - 2 - a n b n -i b n S n -2, (5) 
PiPoPs = a n + b n -i b n S n -i - a n b n -i b n S n - h (6) 
PiPoP2P3 = a n b n + a n b n -i + b n -i b n - 2a n b n -i b n . (7) 

Writing the two-terminal reliability Rel^ 5 (So — > S n ) = S n S n , we immedi- 
ately see that 

S n = Pi Po S n -i(a n -i -> a n -i S n -i p 2 , b n -i -> p 3 ), (8) 



because the transformed graph remains essentially a Brecht-Colbourn ladder, 
provided that some edges are renormalized (fo n _i must be replaced by p 3 , S n -i 
by po, and a n -i by a n -i S n _ip 2 ). 

In the case of imperfect nodes, the assumption 

S n S n -1 + Pn b n S n _i + 7 

n Q"n "n S n -1 (9) 

gives four parameters, for which we want to find a recursion relation. When 
we use eqs. (J5])-([7j) in eq. (jUJ), there is one slightly tricky point: we are actually 
dealing with Boolean functions, so that S^_ 2 = Sn-2 (see for instance [3]) and 

Po a n -i S n -2 (a n -i S n -i P2) 

= Oln-l S n -l S n -2 O-n-1 ipn + a n b n -l S n -2 ~ 0-n b n -l b n S n -2) 
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where the 4x4 transfer matrix M n _i is given by 



M n -1 









S n -2 


\ 


Q-n-1 b n -l S n -2 


a n -i b n -i S n -2 





Q>n-1 b n -l S n -2 


0-n-l 


a n -i S n -2 


b n -l S n -2 


a-n-i b n -i S n -2 


—Cbn-1 b n -l S n -2 


—a n -i b n -i S n -2 


—b n -l S n -2 


Q>n-1 (1 — 2 b n -l) S n -2 ) 



.(12) 



For n = 2, the two-terminal reliability is easily solved because we have series- 
parallel graph: 



Rel^^So — * S 2 ) = S 2 (a 2 + h S 1 b 2 - a 2 h S 1 b 2 ) S , 



(13) 



which leads to a 2 = So, a 2 — 0, j3 2 — b\ So, and 72 = —b\ So, or equivalently 
to cki = 5i = 71 = et (3\ — 1 (see eq. fTT2l ). The final result is 



Relf c) (5 -> S n ) = S n (0 1 0) ■ M n ■ M n _ x 



Mi 





1 

V / 



(14) 



Note that in M\, a\ is quite arbitrary, so that it can be set equal to zero 
without loss of generality. Equation ffT4|) also holds for n equal to 1 or even 
(following a frequent convention that a product of zero matrices is the identity 
matrix); it has also been independently checked for the first values of n, using 
a sum of disjoint products procedure. 

In the case of perfect nodes, the ansatz should be S n = a n a n + (3 n b n + 1n «n b n - 
The relevant 3x3 transfer matrix M'„ is now 



M' r 



Q"n bn 


1 


b n 


n n 


b n 


&n bn 




~b n 




(l-2 6„) 



(15) 



from which we deduce 



Relf C) (,So -> S n , perfect nodes) = (010)- M' n ■ M' n _ x • • • M\ ■ 



v / 



.(16) 



S 



3.2 Generalized fan 




Fig. 5. Triangle-star transformation for the generalized fan. 

Let us now turn to the second architecture. The triangle-star transformation, 
applied to the generalized fan, is displayed on Fig. [51 and we expect the com- 
mon node T to play a special role in the recursion relation. We can write 
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Relj; (So — > S n ) = S n <S n , and see immediately that 

S n = PlP0«S n _l(a„_l -> S n _ip 2 A-l ~> P3), 



(17) 



with 



Pi Po P2 = a n + K-i b n T -a n b n -i b n T, 
Pi Po P3 = b n + a n b n -i S n _i - a n b n -i b n S re _i, 
Pi Po P2 P3 = a n b n + a n b n -i + 6 n - 2 a n 6„ 

Our ansatz is now 

S n = S n a n a n + a n a n T + (3 n b n T + 7 n a n 6 n T. 



(18) 
(19) 
(20) 



(21) 



Here again, we must be careful with the Boolean variable T. When using 
eq. (118j) in eq. (12 II) . all occurrences of T 2 must be replaced by T . The recursion 
relations take the now familiar form 



ft, 

ft 
\7n/ 





-A 


a n - 


i 


Pn- 


-l 







(22) 



where the transfer matrix M n is 
/ 



Mn 



Q"n S n 
















b n S n 


b n S n 


G^n b n S n 


(^n bn S n 


1 


(^n bn S n 


^ &n bn S n 


Q"n bn S n 


b n S n 


a n (1 - 2b n ) S n 



(23) 



When n — 1, we have a series-parallel graph: 

Re4 fan) (S -> Si) = Si(ai + b b x T - a x b h T) S , 



(24) 
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from which we deduce the final expression 



Rel? 11 ^ -> S n ) = (1 T 0) • M n ■ M n _ x ■ ■ ■ M 1 ■ M • 








(25) 



with the additional convention a = 1 for M . We have also checked the 
correctness of eq. (1231) by a sum of disjoint products procedure for the first 
values of n. It also agrees with the expression of Aggarwal et al. [3j for n — 3. 

In the case of perfect nodes, which is often considered in the literature, the 
ansatz S n = a n a n + (3 n b n + j n a n b n leads to 



RelJ^So -> S n - perfect nodes) = (10 0) • M' n ■ M' n _x ■ • • M~\ ■ M' c 



(26) 



V / 



with 





/ 

On 






(In bn 








ct n b n 


1 




(27) 






-6n 


An 


(i-2 6„; 







and again ao = 1 for M . 

It is worth noting that while the sizes of M n and M n for the Brecht-Colbourn 
and fan cases are identical, their matrix elements are clearly distinct. This 
actually leads to strong differences, as will be shown in the following section. 



4 Identical reliabilities p and p 



4-1 Introduction 



While eqs. (fT2l) and (fT4j) give the two-terminal reliability for the general 
Brecht-Colbourn ladder, and eqs. (1231) and (1251) for the fan, we consider here 
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the special case where all edges and nodes have reliabilities p and p, respec- 
tively. Because identical edge reliabilities are usually taken for granted, our 
exact results may help to demonstrate once again the underlying connection 
between combinatorics and reliability theory, most particularly in the enu- 
meration of self-avoiding walks [26f 351152) . This assumption greatly simplifies 
the problem, since all transfer matrices, with the exception of Mq in the case 
of the generalized fan, are identical to a matrix M(p,p) (or M(p,p)): the 
two-terminal reliability is essentially the n th power of one matrix. While the 
determination of its eigenvalues and eigenvectors is one way to solve the prob- 
lem, we adopt in the following the generating function formalism, which is a 
fundamental tool in combinatorics [65j, and makes for a very compact synthe- 
sis of the results. The generating function G(z) is defined by 



(JXJ 

G(z) = ^^t\v,p) z n - 



(28) 



n=0 



From the finite-order recursion relations obeyed by Rel^, the generating func- 
tion is obviously a rational function of z, namely G(z) = N(z)/D(z). In the 
present study, its denominator is at most of degree 4 in z, because in the 

£ 4 -P car ac(lA), where P carac is the 



absence of further simplification, D(z) 



characteristic polynomial of the transfer matrix. The determination of N(z) 
is then very straightforward: we only need to multiply the first terms of the 
expansion of G(z) by D(z) to see the numerator N(z) emerge. The partial 
fraction decomposition of G(z) finally leads to a compact, sometimes very 
simple, analytical expression of RelJ i(p, p), for arbitrary n. 



4-2 Brecht-Colbourn ladder 



4-2.1 Imperfect nodes and edges 

When edge and node reliabilities are p et p, the transfer matrix of eq. (I12p is 
equal to 



M(p,p) 









p 


\ 


P 2 P 


P 2 P 





v 2 P 


p 


VP 


VP 


v 2 P 


-P 2 P 


2 

-VP 


-VP 


p{\ -2p)p) 



(29) 
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and 



Re\ 2 (S ^S n ) =P (0010)- M(p,p) n 





1 

V / 



(30) 



The characteristic polynomial of this matrix is 



carac [P, P\ X ) = X ~ P (2 ~ p) p X - p p (1 - p (p + p - p )) X 

+ (1 - p) (1 - pp) p 2 p 2 x - (1 - p) (1 - p) p 4 p 3 . (31) 

The eigenvalues of M(p, p) are given by the roots of this polynomial of degree 4. 
From the first few values of Rel 2 (S , — »■ S n ) (taking n = 2, 10, for instance), 
we can easily deduce the generating function G#q(z). 



2 2 ly BC\ Z ) 



N, 



Dbc(z)' 



(32) 



with 



4%) 



1 + p (1 — p) p + p 2 p(l — p p (2 — p)) z — p 2 p 2 (1 — p) 2 z 2 
+p 4 (l-p)p 3 (l-p)z 3 (33) 
1— p(2 — p) pz — pp(l — p(p + p 2 — p 3 )) z 2 
+ (1 -p)(l-pp) p 2 p 2 z 3 -(l- p) (1 - p) p 4 p 3 z 4 (34) 



Adding p + p p 2 z (basically, the two-terminal reliability of a single node, or 
that of two connected nodes) to G^q(z), the numerator further simplifies, so 
that we may use another generating function, which of course provides the 
same coefficient of z n for n > 2, namely 



p (1 — p) p z + p 3 (1 — p) p 2 z 2 
DbW) 



Gbc{z) = p 



(35) 



Nbc( z ) wm be the numerator of Gbc(^) i n e Q- d35D- 

What is the nature of the roots ? Numerically, two situations occur when 
< p < 1 and < p < 1: there are either (i) four real roots or (ii) two 
real roots and two complex roots. For instance, when p = 0.9, there are four 
real roots if p < 0.5533938..., but only two otherwise. The separation of the 
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two domains occurs when one root is degenerate. When this happens, the 

dP( BC Hn p- x) 

equalities P^(p,p;x) = and — ca ^ 1 = lead to a polynomial 

constraint satisfied by p and p, given in the appendix (eq. (1A.5j) ). For a given 
p, there exists a unique solution in the range < p < 1, p C rit(p), which is 
displayed in Fig. El 



P. 




4 real roots 



0.2 



m 

P 



0.Q 



o.a 



Fig. 6. Separation of the domains (p,p) for which four real roots exist (under the 
curve) or two real and two complex roots coexist (above it). 

From the inspection of eq. ( |35l) . we may wonder when numerator and denom- 
inator of the generating function may share a common root. A Grobner basis 
calculation shows that this happens only when p = 2, outside our range of 
interest. Therefore, there are four distinct roots most of the times, and when 
there is a double (real) root of V(p C rit(p), P', x), it is indeed degenerate (there 
is no simplification of Gbc(^))- The exact expression for Rel2(So — > S n ) can 
then be deduced using a partial fraction decomposition of Gbc( z )- When all 
the eigenvalues Aj are distinct, 



Gbc(z) 



E 

A, 



1 



T- 1-XiZ' 



(36) 



consequently, the two-terminal reliability reads 
Re\ 2 (S ^S n )=^2 on K- 

A* 



(37) 



In the Brecht-Colbourn case, the root A max of greatest modulus is always real. 
Unless p and p are equal to 1, A max < 1, so that Rel^^So — * S n ) decreases 



14 



as n 



oo, essentially as a power- law behavior 



Rel 



(BC), 



*2 '(S — > S n ) ~ A" ax . 
4-2.2 Imperfect edges and perfect nodes 



(3* 



We devote this section to the case of imperfect edges and perfect nodes, since 
the Brecht-Colbourn ladder was initially investigated in this configuration 
|13|[T4"] . We immediately find 





f p 2 


1 


p 2 ) 


n 


(0) 




Ref c ' p=1 \S - S n ) = (0 1 0) ■ 


p 


V 


p 2 




1 


, (39) 




K -p 2 


-p 


p(l-2p)j 









so that the 25-node case HU leads to 



Rel 2 (S" -» 5 24 ) =P (1 + 78p + 1121/ + 6633 p 6 + 11554 p 4 - 57525/ 

-279450/ + 89578 p 7 + 2570040 p 8 + 1431183 p 9 



-17159566p 
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12166498 p 11 + 98985590p 12 
+33119917p 13 - 495566666 p 14 + 212008622 p 15 
+1867178285p 16 - 2888906214 p 17 - 3066846055p 



18 



+14427083319p 19 - 14178781875 p zu - 17955754991^ 



20 



,21 



+80808979717p 



22 



144754404751 p 23 + 174732303288 p 24 



-158925117297p 25 + 113702258108 p 2e 



65190712312p 



27 



+30133254848p 



2N 



11 197537798 p 29 + 3308571601 p 30 



761406139 p 31 + 131805146 p dZ - 16170009p 



32 



v 33 



+1254886 p 34 - 46368 p 



35 ^ 



(40) 



The coefficients are already quite large, even though the number of nodes 
remains limited. To the best of our knowledge, the exact two-terminal reliabil- 
ity polynomial was calculated only for a 10-node ladder [13J; we recover this 
result by taking n = 9 in eq. ( 1391) . 

The three distinct real eigenvalues A» (i = 1,2,3) of the transfer matrix ap- 
pearing in eq. ( |39l) . one of which is negative, are displayed as a function of p 
in Fig. [71 



Their analytical expression is (£ E {1, e 2t7T ^ 3 , e 2i7r / 3 }) 



\ = ~p(2-p) + - 



'A + IV27VB 



1/3 



1/3 



(41) 



15 



0.4 0$ 

p 



Fig. 7. Variation with p of the three eigenvalues for the Brecht-Colbourn ladder 
with perfect nodes (see eq. {H}). 

where 



A=-p 2 (9 - 43 p + 60 p 2 - 39p 3 + lip 4 ) , (42) 
£= (1 - p) 2 p 3 (4 + 9p + 16p 2 - 88 p 3 + 98 p 4 - 32 p 5 - 8p 6 + hp 7 ) ;(43) 

note that i3 is always positive for < p < 1. 

The generating function for the Brecht-Colbourn ladder with perfect nodes is 
simply 



4 p c =1) (*) 



1 — p (1 — p) + p 3 (1 — p) z A 



1 — p (2 — p) z — p (1 — p) 2 (1 + p) z 2 + p 2 (1 — p) 2 



(44) 



From its partial fraction decomposition, which is necessarily of the form G# c ^ (z 



8=1 - 1 - 



A,; 2 



because the eigenvalues are distinct, we deduce after simplification 



Rel 2 (S ^S n ,)=E 



p(l-p) 2 -(l-p) \ t -Xj 



\ 3p (1 - p) 2 - 2(1 - p) 2 (1 + p) Aj — (2 - p) A? 



\».(45) 



From eq. (14"5]) . it is clear that the large n limit of Rel 2 (S , — > S n ) is a power- law 
behavior, in which the eigenvalue of maximum modulus quickly prevails over 
the others when p is close to unity, even for moderate values of n. 
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Table 1 

Values of the F^s for the two-terminal reliability polynomial of the 25-node Brecht- 
Colbourn ladder. 



i 


77 

Fi 


i 


77 

Fi 


i 


771 


i 


771 

Fi 





1 


9 


1125395882 
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762855455898 
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12199394435 


20 


725278875430 
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16103 
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o o 'm oor /i /i o t 

32708854487 


21 


562806091836 
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30 


13648753 
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175418 
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76833130394 
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371300292894 


31 


1334810 
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1484837 


14 


158368734141 


23 


206539411448 


32 


93929 


6 


10151340 


15 


286502593795 


24 


96061397122 


33 


4368 
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57524387 


16 


454444238576 


25 


37052347922 


34 


113 


8 


275139029 


17 


630595957484 


26 


11756780232 


35 


1 



2n-l 

We can also write Rel 2 (5o -> #„) = p 2 71-1-1 (1 - p)\ a well-known 

i=0 

expansion of reliability polynomials [26]. The i^'s for the 25-node ladder are 
given in Table [1] (the maximum value was equal to 8078 in the 10-node ladder 
|13j). Let us call m = n + 1 the number of nodes. Using linear regressions on 
the first values of Re^So — > S n ), it is straightforward to find 



Rel 2 m odd \p)=p"' 2 ^ (1 — p)' 2 " H p 2 

m 4 + 72 m 3 + 350 m 2 - 2376 m + 2337 m+ 3 . 
H — P 2 (1— p) 



m 2 - 12 m - 21 



;i-p) : 



384 



(46) 



t-j .(mevcn), s Tt m 3m^6 (m - 2) (m 2 + 38 771 + 24) m+2 

Rel 2 \p) = — p 2 (1 -p) 2 + - ^~~48 -P 2 (1-p) 



m-2 m 4 + 122m 3 + 2304 m 2 - 5472m- 13440 m+4 , 

1 12 -» 2 (1-p)" 

3840 F V ^ ; 



+ 



Finally, the comparison of the exact results with the various lower bounds 
proposed in [2] is given in Table [2] and Fig. El It shows that the lower bound 
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Table 2 

Comparison of various lower bounds with the exact two-terminal reliability of a 
25-node Brecht-Colbourn ladder [14]. 



V 


Kruskal-Katona 


MinCost(edp) 


Br echt- Colb our n 


exact 


0.75 


0.031682 


0.054681 


0.054681 


0.625163 


0.80 


0.068803 


0.119917 


0.430912 


0.773696 


0.82 


0.092654 


0.161200 


0.558991 


0.824038 


0.84 


0.124041 


0.214282 


0.669269 


0.867950 


0.86 


0.165305 


0.281396 


0.761945 


0.905042 


0.88 


0.219694 


0.364529 


0.837486 


0.935251 


0.90 


0.291856 


0.464826 


0.896659 


0.958806 


0.91 


0.336579 


0.521297 


0.920440 


0.968231 


0.92 


0.388392 


0.581555 


0.940574 


0.976194 


0.93 


0.448415 


0.644934 


0.957259 


0.982785 


0.94 


0.517724 


0.710375 


0.970720 


0.988109 


0.95 


0.597041 


0.776313 


0.981207 


0.992275 


0.96 


0.686113 


0.840514 


0.989003 


0.995400 


0.97 


0.782518 


0.899895 


0.994420 


0.997608 


0.98 


0.879474 


0.950274 


0.997801 


0.999024 


0.99 


0.961964 


0.986085 


0.999524 


0.999778 



of Brecht and Colbourn is rather good for p close to unity, but its sharpness 
decreases for p < 0.8. 



4-2.3 Imperfect nodes and perfect edges 

The configuration of imperfect nodes and perfect edges, although less stud- 
ied than the previous one, has nonetheless been considered in several pa- 
pers [Tlll2ll35j . The Brecht-Colbourn ladder case has been studied in detail by 
Graver and Sobel [35]. The relevant transfer matrix is now equal to 
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M(p=l,p) 



(48) 
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Fig. 8. Comparison of the exact value (^^— ) with lower bounds given in ref. 
for the 25-node ladder:(A) Kruskal-Katona (■) Min Cost (♦) Brecht-Colbourn 



Its characteristic polynomial is x 2 (x 2 — px — p(l — p)), which simply leads 

to the eigenvalues — which plays no role in the final expression of the two- 

1 



terminal reliability — and \± = — (p±-y/4p — 3p 2 ). The generating function 
is found to be 



1 — p z — p (1 — p) z 2 A+ — A_ \1 — A + 2 1 — X- z / 



,(49) 



which gives 



Rel 2 (5'o -> S n 



V4: P -3p 2 



(a 



n+l \ n +i 



(50) 



This is almost Graver and Sobel's result [35], obtained through a combinatorial 
argument. Their final expression differs from ours because their source and 
destination (in our notation, So and S n ) are perfect. 
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4-3 Generalized fan 



When the reliabilities are p et p, the transfer matrix of eq. (1231) is equal to 



M(p, p) 



VP 



p 2 p 



\-p P 



VP 



v 2 P 



2 

-VP 



VP 



-VP 



v 2 P 



v 2 P 



\ 



p{l -2p)p) 



(51) 



and 



Rel 2 (S -> S n ) 



-(lpOO) • M(p,p) n+1 
V 








(52) 



the 1/p prefactor being a consequence of the condition a 
characteristic polynomial of this matrix factorizes nicely: 



± carac IF ' r ) 



■T) 



(x — 1) (x — p p) [x — p (1 — p) p 



1 in M . The 



(53) 



so that the eigenvalues are 1, pp and p{l — p) p, the latter being of degree 
2. The presence of 1 among the roots should not be surprising. Indeed, even 
when n goes to infinity, the two-terminal reliability between So and S n is larger 
than pp 2 , in stark contrast to the Brecht-Colbourn case, where the reliability 
vanishes. Since a power-law behavior is still expected, the only possibility is 
that the largest eigenvalue is of modulus 1. 

The generating function G^ n \z) is derived using the recipe described in the 
Brecht-Colbourn case. The additional convention Rel2(So — > So) = p leads to 



G (fan) (^ 



N^ an \z) 



p) (2 - pp)]+ppz 2 [(2 + pp){l-p) 



£>(fen)(*)' 

1 - z [1 +pp{l 

+V 2 P(V~ P)] ~ V 2 (1 ~ V) 2 P 2 z 3 , 
D^ Q \z) = (l -z) (1-ppz) (l-p(l-p)pz) 2 . 



(54) 

(55) 
(56) 
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Because all the roots have simple expressions, we expect to find a simple, 
analytical expression for Rel^ (5o ~^ S n ). The partial fraction decomposition 
of eq. (153|) gives indeed 



G fr*)( z ) = P 2 (1~P) l-2p(l-p)p + p 2 p 2 (l-3p + p 2 



l-p(l-p)p) l-p(l-p)pz 
2 l-pp(2-p) 1 



+pp 



l-p(l-p)p (l-p(l-p)pz) 
p 2 p 3 1 . 



+ p(l-p) . (57) 

(l-p(l-p)p) 2 1-z ^ V 1 1 l-ppz K ' 

The new feature of eq. (jSTj) is the (1 — p (1 — p) pz)~ 2 term. Because 

1 



z yi = E("+ 1 )* B * B > 

1 n=0 



the final expression of the two-terminal reliability is 



Ret n \So^S n ) = +p n (l-p) n p 



n _n+2 



np 



pp(2 -p) 



l-p(l-p)p 

1 — p p (2 — p) + p 2 (1 — p) 2 p 2 " 



p{l-p) p)' 



+p n p n+1 (1 - p) + 



p 2 p 3 



;i -p(i -p) p) 



2 ' 



(59) 



where n appears in a prefactor, not only as an exponent. The second term of 
eq. (159]) is the asymptotic limit when n — > oo. We see that the reliability of 
the path #0 T — > S n is enhanced by 1/(1 — p (1 — p) p) 2 . When nodes are 
perfect, p must be set to one in eq. ( 1591) . and the contribution of the eigenvalue 
pp vanishes (the transfer matrix of eq. (|2"7|) is 3 x 3). This leads to 



Rel (fan;p=l) 



(S ^S n )=p n (l-p) 



n+2 



np 



+ 



;i + p 2 



i-p(i-p) (i-p(i-p)y 



p 



p(i -p)) 



2 ■ 



(60) 
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5 Zeros of the two-terminal reliability polynomials 

The structure of the different reliability polynomials may be understood by 
studying the locations of their zeros in the complex plane. Such a study has 
been fruitfully performed in the case of the chromatic polynomial [91lll|l57j . 
most notably in the context of the four-color theorem. In reliability studies, 
some effort has been done to discover general properties for the all-terminal 
reliability Kel^ip) [2T|28ll51j . its main byproduct being the Brown-Colbourn 
conjecture [16], according to which all the zeros are to be found in the re- 
gion |1 — p\ < 1. Although valid for series-parallel graphs, this remarkable 
conjecture does not strictly hold in the general case (but not by far) [55] . As 
mentioned in the introduction, the all-reliability polynomial is linked to the 
Tutte polynomial, an invariant of the graph. It has also been studied exten- 
sively by Chang and Shrock for various recursive families of graphs [19], who 
give the limiting curves where all zeros of the polynomials converge. 

In this section, after briefly recalling general results of the literature, we study 
the roots of Rel 2 (jo, p) = in the complex plane for a fixed p, in order to 
see whether some insight may also be gained in this case. Admittedly, this 
polynomial depends on the couple (source, terminal), but structures are still 
expected. We show that the zeros tend to aggregate along portions of alge- 
braic curves, which can substantially differ even for the two families of graphs 
under consideration, even though they have the same all-terminal reliability. 
Moreover, structural transitions occur at p = \. 

5.1 Calculation of the limiting curves (generalities) 

As n grows, the number of zeros of the reliability polynomial in the complex 
plane increases. Because of the matrix transfer property, we have recursion 
relations between relability polynomials corresponding to successive values of 
n. The general treatment of the problem has been done by Beraha, Kahane, 
and Weiss [H], but may be understood in the following, simplifying way: if the 
reliability polynomial is of the form J2% a i \ {p) n (where Aj are the eigenvalues 
of the transfer matrix), then at large n, only the two eigenvalues of greater 
modules, say Ai and A 2 , will prevail, so that the reliability polynomial will 
vanish when |Ai(p)| = | A2 (j>) | (of course, it might be three or more eigenvalues 
of equal modulus; the present oversimplification works quite well here). This 
defines a set of curves in the complex plane, where all zeros should accumu- 
late in the n — > 00 limit. The interested reader should refer to the work of 
Salas and Sokal [57] for a very detailed discussion of the convergence to the 
limiting curves. This behavior is not modified when one or more of the a^s is 
a polynomial in n [S], as will become apparent for the generalized fan. 



22 




Fig. 9. Localization of the complex roots p of the two-terminal reliability polynomial 
for a Brecht-Colbourn ladder of 151 perfect nodes and 299 edges. 

5.2 Calculation for the Brecht-Colbourn ladder 

The recursion relation obeyed by Rel^ 5 (So — * S n ) = Rel^ is easily deduced 
from eq. ( 1311) ; it reads 

Rel^ip) =p(2-p)p R4 n_1) + pp{l-pp-p 2 p + p 3 p) Rel^ n ~ 2) 

-p 2 (1 - p) p 2 (1 - p) Relf" 3) + p A (1 - p) p 3 (1 - p) Rel^ n - 4) . (61) 

and leads to very quick calculations using mathematical softwares such as 
Mathematica [15] . Note that this recursion relation of order four becomes of 
order three when p — 1. We can then look for solutions of Rel^Qo) = in 
the complex plane. Of course, the degree of the polynomial, as well as the 
magnitude of its coefficients, will increase with n. 

5.2.1 Perfect nodes 

Let us begin by setting p = 1. A first step is to display the roots of ReL^Qs) = 
0. These zeros have been calculated using the NSolve [] routine of Mathe- 
matica; because the polynomial coefficients can be very large, a numerical 
accuracy of several hundred digits is sometimes necessary. Figure [9] shows the 
location of the zeros of Rel^ 1 (p) in the complex plane (it corresponds to a 
graph with 151 nodes and 299 edges). 

On the right half of the complex plane, we have a simple, open curve crossing 
the real axis at D and extremities at C and C*. On the left half, however, we 
have what looks like a closed curve — intersecting the negative real axis in A 
and B — as well as some dots on the negative real axis, between A and the 
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origin. 



As mentioned above, it is well known that as n increases, the zeros accumulate 
at particular locations, constituted by segments/portions of algebraic curves. 
For n = 150, these limits are almost reached, even though the "sampling" is 
not uniform. We can calculate the limiting curves by considering that the roots 
of the characteristic polynomial are \e id/2 , \e- i9/2 and /i. The coefficients 
of the characteristic polynomial imply constraints between A (which may be 
complex), p and 9, or more precisely, t = cos | or T = cos 9, which will be 
used in the following. For instance, p = (2 — p)p — 2 A cos |. After eliminating 
A and p, we obtain a compatibility condition which must be satisfied by p, 
and T (the polynomial 7-3 (p, p, T) is given in eq. fl A. 3[) ) : 



P 3 (p, p = 1,T) = 2 + 9p - Ap 2 - 30 p 3 + 38 p 4 - 10p 5 - Qp 6 + 3p 7 

+ (2 + Up - 8p 2 - 34 p 3 + 48 p 4 - 18p 5 - 2p 6 + 2p 7 ) T 

-Ap(l-p) 2 (l-p + p 2 ) T 2 — 8p (1 — p) 2 T 3 
= 0. (62) 

The additional requirement \p\ < |A| translates into an additional contraint 
between p and T. 

A plot of this parametric set of curves shows indeed that the whole portion of 
the negative real axis between A and the origin is indeed a solution, so that 
we could expect more zeros there when n > 150. 

A few critical points can also be deduced from eq. (162]) . For instance, D cor- 
responds to a real solution of Vs(p, p = 1,T = —1) = 0, namely pr, = 1+ 2 V ^ . 
C and C* are points of the complex plane defined by pc and its complex 
conjugate, which are roots of P3Q9, p = 1, T = +1) = 0, that is 

A + 9p + 16p 2 - 88p 3 + 98p A - 32p 5 - 8p 6 + 5p 7 = 0. (63) 



This polynomial already appeared in the expression of B. Note that there 
are more than one pair of complex solutions to this equation, the relevant 
one is given by pc ~ 1.011578 + i 0.607394. The determination of the com- 
plex numbers pa and ps associated with A and B is a little more elaborate. 
Both pa and ps are real and negative; they correspond to roots of degree two 
of eq. fl62l) . In order to find them, we must have V3 (p, p = 1 , T) = and 

(p, p , — ) _ q Thggg t wo relations are satisfied for special values of p 

op 

(or T). After substitution and elimination of T using Mathematica, we obtain 
a product of polynomials, which must cancel for p = Pa an d p = Pb- For the 
sake of completeness, we provide the needed polynomials in appendix [A] A 
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close numerical study allows to identify the relevant polynomials: 

• p A « -0.2879878 is solution oiVi(p,p= 1) = (see eq. (TO]) ), that is 

= 1 - p- lip 2 + 15p 3 - 3p 4 - 2p 5 -p 6 +p 7 , (64) 
attained for T « 0.138176, a root of 

= -1 + 11T-27T 2 -26T 3 + 140T 4 + 240T 5 + 144T 6 + 32T 7 .(65) 

• p B « -0.1849482 is solution of V 2 (p, p = 1) = (see eq. (TAT2]) ). i.e., 

= 1 - 4p - 9p 2 + 101 p 3 - 413 p 4 + 1019 p 5 - 1761 p 6 

+2151/ - 1864 p 8 + 1097/ - 386 p 10 + 60 p 11 , (66) 

attained for T m —0.9511957, one of the solutions of 

= 3040707 + 12576464 T + 15322821 T 2 - 4376828 T 3 

-22559186 T 4 - 101 12842 T 5 + 9510320 T 6 + 7762048 T 7 
-1337920 T 8 - 2068608 T 9 + 8192 T 10 + 204800 T 11 (67) 

5.2.2 Imperfect nodes 

We can perform the study for p < 1 by following the same lines as above. 
Although more cumbersome because there are now four possible roots instead 
of three, the calculations are interesting nonetheless because we may vary one 
parameter, p. A natural question is: can it affect the global structure of the 
zeros ? The answer is yes, and the location of the zeros undergoes a kind of 
"structural transition" for a particular value of p, namely |. 

Let us first consider p > |. The location of the zeros is then qualitatively 
similar to that described in Fig. [91 the whole structure merely expanding 
from the origin. After comparison with the nu merica l values for n = 150, 
D appears to be associated with po = | (1 + yT + the relevant solution 
of Vzip-, p, T = —1) = 0. Likewise, C and C* are given by pc, which is one 
complex solution of Vs(p, p, T = +1) =0 (V3 is given in eq. ( 1A.3I) ). 

The determination of pa and pb is more tedious, because Vs{p, p, T) is now a 
polynomial of degree 6 in T. Here again, we must find the common zeros of 
Vs(p, p, T) and its derivative with respect to p. The elimination of T, performed 
using Mathematica, leads to a polynomial in p and p, which must vanish. 
Actually, this polynomial can be factored, and a numerical comparison with 
the zeros obtained for n = 150 shows that p A is solution of V\(p, p) = (see 
eq. (lA.ip ). Similarly, p B is a solution of ^(p, p) = 0, where P 2 is given in 
eq. d£2}. 
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Fig. 10. Localization of the complex roots p of the two-terminal reliability polyno- 
mial for a Brecht-Colbourn ladder of 151 nodes and 299 edges, when p = 10~ 2 . 

For p = i, p A « -0.4359355, p B « -0.2885759, p c « 0.748541 + i 1.03759, 
and = 2. Note that and ps are zeros of polynomials in p of degrees 
equal to 22 and 30, respectively (having p ^ 1 does not simplify the problem). 

For p <\, the above expressions for p A , Pb, and pc are still valid, even though 
C and C* may now belong to the left half-plane. However, on the right half- 
plane, the structure of zeros undergoes a drastic transformation, as shown in 
Fig.[10]for p = 10~ 2 . D is now an angular point, and zeros may be found on the 
real axis between D and E. After numerical comparisons between structures 
of zeros and limiting curves have been performed, it appears that pn is now 
a real root of Vi(p, p) = (actually, the third one, in decreasing order), while 
Pe is another real root of V%{p, p, T = +1) =0. 



5.2.3 Asymptotic limits when p — > 

In this section, we address the asymptotic dependence as p decreases to zero 
of all the critical points A, B, C, D, and E. A first step consists in finding 
the corresponding p's for (very) small values of p. It is not too difficult to 
observe that these numerical values have a p -1 / 3 dependence. We can then 
use improve on this knowledge by using the algebraic equations of which they 
are solutions, making again use of Mathematica for the asymptotic expansions. 
A comparison with the numerical results determines the true leading term of 
the expansion; after some work, we obtain 




1/3 



35 + 12 y/E 
90 



+ 




(68) 
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Pb(p)=-(jJ ' + « + 0(p 1 / 3 ), (69) 
where 

1 f/2531 + 15v / 3l593\ 1/3 /-2531 + 15 v / 31593\ 1/3 

x = i6 U 2 j i 2 j " 8 < 70 » 



is the real root of 5 x 3 + 8 x 2 + 8 x — 1 = 0. Numerically, x ~ 0.11166155366 
and k = x 1 / 3 ~ 0.48154242495. The next-order calculation provides 



K 1 70147451 k 2 + 22890531 k - 3689542 1/3 
PB W-~pm + 50 1685743 + 778683 k- 121256 (P ^ ( ^ 

S v ' 

a 

0.48154242495 



pi/3 



+ 0.38969988720 + 0(p 1/3 ). (72) 



Proceeding in a similar way for the critical points defining the rightmost struc- 
ture in the complex plane, we find 



pM= (l^r + }^ + 0^ ) (73) 



2p J 30 

/3-v/sV 73 2 /3 + v/oY 76 
Mp) = — = +77 S V -75 + 35 V5 



2p J 15 V 2 / v pV6 
+i^ + 0(^) (74) 

V3 n ,„ , /fx 5/6 , ei7r/3 



2p / 15 V 2 / v pV6 
0(p 1/6 ) (75) 
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Note that the expansion pc{p) is nothing but Pe{p), with the transformation 
p — ► e _2l7r p. It corresponds to the second complex root, by decreasing order 
of the real part. 



5. 3 Calculation for the fan 



We expect here much easier calculations, since all the eigenvalues are extremely 
simple: 1, p p, and p (1 — p) p. They are, indeed, since the limiting curves can 
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be expressed analytically (see Table [3]). It is worth noting that nonetheless, 
other structure transitions occur at p = 1 and p — ~. Using the relevant 
recursion relation between successive reliability polynomials, we can again 
calculate Rel^^So — > S n )(p) rapidly using Mathematica. 




-n.s 



Re(p) 

Fig. 11. Location of the roots of Rel^™ -150 ^ (p) = for the fan, with p = 1. 



For p = 1, we see on Figure [TT] that the limiting curve is now essentially a 
closed one, even though the zeros appear to belong to two different sets which 
slowly join as n increases. The presence of n in a prefactor as in eqs. (1591) 
and (l60j) does not affect the general definition of the limiting curves [8]. The 
limiting curve is very easy to obtain, since we have only two eigenvalues, 1 
and p (1 — p). For them to have the same modulus implies p (1 — p) = —e 10 , 
so that we get p = \ (l ± \f\ +4e if? ) . 

When p is strictly less than 1, the whole structure is altered, as witnessed by 
Figure d2] for p = 0.9999. Portions of circles centered at the origin and at the 
point (1,0) are added. The reason for such additional structures comes from 
the existence of a new eigenvalue pp. They persist down to p — \. For p < \, 
the only eigenvalues that matter are actually 1 and p (1 — p) p, so that the 
structure looks again very much like Fig. [TH 



Their derivation being straightforward, we simply give the analytical expres- 
sions of the limiting curves in Table [3j When p — > 0, the structure expands 
with a scaling factor of p -1 ^ 2 , in contrast with the p" 1 / 3 obtained for the 
Brecht-Colbourn ladder. 
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Fig. 12. Location of the roots of Rel^T (p) = for the fan, with p = 0.9999. 
Table 3 

Analytical expressions of the limiting aggregation curves for the location of p such 
thatRe4 fan) '™(p,p) = 0, as a function of p, the node reliability. 



node reliability 


parametrization 


validity range 


p=l 


p= \ (l± v / l + 4e i6 ') 


-1 < cos 6 < 1 






-1 < cos 9 < 1 


\ <p< 1 


p=!(l + ^l + i^) 


-l<cos#<^ (^-3) 






cos 9 > 7T- 

z p 




p= l + e ie 


COS > — 1 


< p < \ 




-1 < COS# < 1 



6 All-terminal reliability 

As mentioned in the introduction, the all-terminal reliability Rel^ is another 
useful measure of the network availability, giving the probability that all nodes 
are connected. In this case however, the node reliabilities are a mere overall 
factor, so that they may be considered equal to 1 for all practical purposes 
[26] . Chang and Shrock [19] gave the explicit expressions of ReU for various 
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recursive families of graphs, among which the Brecht-Colbourn ladder, with 
the same reliability p for all edges. The all-terminal reliability of the general- 
ized fan has been calculated by Neufeld and Colbourn [3H] • The results are of 
course identical, because at each step, two edges are added, while the common 
original graph is the triangle (the complete graph K%). For the sake of com- 
pleteness, we slightly generalize their result for edges with distinct reliabilities. 
Here again, the final, analytical expression can be written in a concise form 
using transfer matrices. In the context of graph theory, this can be viewed as 
the factorization of a particular value of the multi-variate Tutte polynomial, 
considered by Wu [73] and Sokal [53]. 




Fig. 13. Deletion-contraction simplification for the all-terminal reliability of the 
Brecht-Colbourn ladder. 



We expect the all-terminal reliability lZ n = ReU(So <-> S n ) for the Brecht- 
Colbourn ladder to exhibit the same behavior as the two-terminal reliability, 
with a transfer matrix depending on a n and b n , and a generating function that 
is a rational fraction of p (the common reliability of links) and z. As mentioned 
above, this calculation should be somewhat easier because all nodes can be 
considered perfect without loss of generality [26]. We can then use the usual 
pivotal decomposition to establish a relationship between lZ n and lZ n -i, as 
represented in Fig. [131 We find 

K n = (l- b n ) dnKn-! + b n TZ n -i(b n -i -> b n -!// a n ) (76) 
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Using the ansatz TZ n = a n (a n + b n ) + 7 n a n b n , we get 



/ 



a. 



( 



i + K-i 


a n-l b n -l 








a n _i (1 - 2 6 n _i) J 




lTn-1 / 



(77) 



with ai = 71 = and a\ = 1. Calling M n _i the transfer matrix of eq. (i77|) . 
we obtain 



ft n = (1, 0) • M n • M n _x • • • Mi 



(78) 



When all reliabilities are equal to p, the transfer matrix is simply 
/, 

M> = 



2p 


^2 


l-3p 


p(l-2p) 



(79) 



P 



The eigenvalues of M' are £± = ^ ( 3 — 2 p ± -y/5 
generating function is 



1 — P (3 — 2 p) z + p 2 ( 1 — p) z 2 



8p + 4p 2 ) , so that the 



50) 



and 



a/5 - 8p + 4p 2 



^•n— 1 £*n— 1 



n > 2 



The last results were obtained in [19], while the same expression can be found 
in the paper by Neufeld and Colbourn |49j . 



7 Conclusion and Perspectives 



We have given the exact solution of the two- and all-terminal reliabilities for 
the Brecht- Colbourn ladder and the generalized fan for arbitrary size and indi- 
vidual element reliability. While simple, these graphs correspond nonetheless 
to realistic network architectures, especially in telecommunication networks 
for IP transport. Node and edge failures are put on an equal footing, and the 
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simple formulae relying on transfer matrices may be directly implemented. We 
have also given the analytical solution of the two- and all-terminal reliabilities, 
when reliabilities are p and p for edges and nodes, along with their rational 
generating functions. The locations of the zeros of the two-terminal reliability 
polynomials differ for the two families of graphs considered in this paper, even 
though their all-terminal reliability is identical. They possess structures in the 
Re(p) < region, and exhibit transitions for particular values of p. 

Even though the delta-star transformation has also been successfully applied 
to the case of simple ladders [BB], it may not be so easy to use in more com- 
plicated networks. However, it seems quite clear that a similar decomposition 
through transfer matrices of two-terminal reliabilities should occur for ladders 
of greater width too. In order to make such calculations useful for applications, 
imperfect nodes as well as imperfect edges must be considered. All is needed is 
a recursion relation between successive graphs, when one "elementary brick" 
is added. This implies a new expression for the deletion-contraction theorem, 
in which the linearity with respect to all individual edge or node reliabilities 
must be preserved, whereas — to our knowledge — edge reliabilities are of- 
ten renormalized to account for the unreliability of the nodes they connect 
[317U|I71] . This new expression will be given elsewhere, along with an applica- 
tion to other recursive families of graphs such as the K± ladder and the K3 
cylinder [57] . 

What should we expect ? Basically, the same kind of behavior as detailed in 
the present work, with a factorization of the reliability in terms of transfer 
matrices, the dimension of which substantially increase to reflect the interplay 
of different edges/nodes in the overall reliability, and the underlying algebraic 
structure of the graph (one cannot escape the intrinsic complexity of the prob- 
lem...). 

There are obviously many directions in which this work may be further ex- 
tended, for instance in the estimates of bounds. The first one is to use the 
present results, or their future extensions [671169] . as possible upper or lower 
bounds to more complex graphs. If a graph under consideration looks like - 
or made as such — a special instance of a recursive family of graphs, we ex- 
pect the generating function to be a rational fraction again. The dimension of 
the corresponding transfer matrix may be probed by trying to find recursion 
relations between successive reliability polynomials. Of course, if the dimen- 
sion of the corresponding transfer matrix is large — which is most likely to 
happen — the degree of the numerator and denominator of the fraction may 
be too large for a complete solution to be obtained easily. Even so, the knowl- 
edge that the true generating function is rational may be an indication that 
an approximate generating function might still be quite useful, because Pade 
approximants [Hl50] are known for their devilish knack of getting very close 
to the exact function. This will be addressed elsewhere [HH]. Our calculations 
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may also provide some information on some combinatorial issues, such as the 
enumeration of self-avoiding walks on lattices of restricted width. 

Other quantities of interest may be deduced from the general expression of 
the network reliability, among them the sensitivity of a equipment [66J, the 
influence of scheduled maintenance on the overall network availability, and 
the failure frequency [371162] , which can all be deduced from various partial 
derivatives of the two- and all-terminal reliabilities [68]. Since the reliability 
of each equipment appears in only one transfer matrix, the computation of all 
these network parameters is straightforward. One should also be aware that 
the reliability of each equipment is not known with absolute accuracy. The 
consequence of this uncertainty on the overall reliability has been addressed 
by Coit and collaborators in the case of series-parallel reducible networks 
[2^125] . Our transfer matrix factorization makes this topic another instance of 
a product of random matrices, definitely a vast field in mathematical physics 

1221 ■ 

Finally, the exact results found for classes of arbitrarily large networks may 
prove useful for testing different algorithms (Monte Carlo, genetic, OBDD, etc) 
in numerically exacting configurations, where edge and node unreliabilities 
must both be taken into account. 
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A Polynomials used in the determination of critical points for the 
zeros of Rel 2 (p) (Brecht-Colbourn ladder) 

A.l Vi(p, P ) 



Vi(p,p)=p-P (l + 9p + 3p 2 ) + p 2 p (41 + 37p- lip 2 ) 
+p 3 p (-39 - 195 p + 99 p 2 - 4 p 3 ) 
+p A p (16 + 147p-411p 2 + 43p 3 + 11 p 4 ) 
-p 5 p 2 (-54 - 2300 p + 123 p 2 + 67 p 3 + 7 p 4 ) 
-p 6 p 2 (126 + 4700 p + 1865 p 2 + 559 p 3 - 59 p 4 + 13 p 5 ) 
+p 7 p 2 (41 + 4859 p + 2324 p 2 + 6400 p 3 + 1135 p 4 + 189 p 5 + 14 p 6 ) 
+p 8 p 3 (-2701 + 4149 p - 11762 p 2 - 8744 p 3 - 2515 p 4 - 311 p 5 + 4p 6 ) 
-p 9 p 3 (-781 + 12546 p - 6054 p 2 - 14564 p 3 - 11548 p 4 - 3012 p 5 - 124 p 6 + 8 p 7 ) 
-p 10 p 3 (96 - 13197 p - 3996 p 2 + 4727 p 3 + 17901 p 4 + 11943p 5 + 1464 p 6 + 4p 7 ) 
+p n p 4 (-7185 - 3814 p - 10626 p 2 + 5695 p 3 + 21337 p 4 + 5871 p 5 + 304 p 6 ) 
-p 12 p 4 (-2040 + 3598 p - 10478 p 2 - 17573 p 3 + 19235 p 4 + 11669 p 5 + 1219 p 6 ) 
+P 13 p 4 (-240 + 6353 p + 3091 p 2 - 30176 p 3 + 6625 p 4 + 14726 p 5 + 2128 p 6 ) 
-p 14 p 5 (3613 + 12361 p - 25852 p 2 - 5221 p 3 + 14188 p 4 + 2147 p 5 ) 
+p 15 p 5 (947 + 11148 p - 14935 p 2 - 8341 p 3 + 9720 p 4 + 2095 p 5 ) 
-p 16 p 5 (96 + 5829 p - 7507 p 2 - 1909 p 3 + 690 p 4 + 3066 p 5 ) 
+2p 17 p 6 (993 - 2513 p + 3568 p 2 - 4056 p 3 + 2046 p 4 ) 
-p 18 p 6 (418 - 3991 p + 10333 p 2 - 10532 p 3 + 3785 p 4 ) 
+p 19 p 6 (41 - 2397 p + 7184 p 2 - 7257 p 3 + 2430 p 4 ) 
+3p 20 (1 - p) p 7 (299 - 684p + 367p 2 ) 
-p 21 (1 - p) p 7 (185 - 519 p + 322 p 2 ) 
+p 22 (1-p) P 7 (l"2p) (16-15p) 

-19 p 23 (1 - p) 2 p 8 + 8 p 24 (1 - p) 2 p 8 - p 25 (1 - pf p 8 (A.l) 



A.2 V 2 (p,p) 



Vz(p, p) = 

9 + 2p (-93 + 37 p) + p 2 (1375 - 1139 p + 235 p 2 ) 
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+p 3 (-4548 + 6619 p - 2306 p 2 + 380 p 3 ) 

+p 4 (7551 - 22762 p + 7011 p 2 - 2415 p 3 + 335 p 4 ) 

+p 5 (-6210 + 53479 p + 2160 p 2 + 7102 p 3 - 1462 p 4 + 154 p 5 ) 

+p 6 (2025 - 85209 p - 79460 p 2 - 18121 p 3 - 11141 p 4 + 223 p 5 + 29 p 6 ) 

*7 /_._. , . O Q . A . . . K\ 



84879 + 277027 p + 106511 /? 2 + 202610 p 3 - 36581 p 4 + 1170 p 5 ) 
-46170 - 521894 p - 587852 p 2 - 1536749 p 3 + 263741 p 4 - 39391 p 5 + 547 p 6 ) 
'-10125 - 611502 p - 1935829 p 2 - 7279427 p 3 + 440373 p 4 



+375501 p° + 2506 p b ) 
+p n p 2 (204255 + 5324974 p + 53128592 p 2 + 28420289 p 3 + 5054943 p 4 
-4185113 p 5 + 12850 p 6 ) 

42525 - 4776213 p - 87000246 p 2 - 104194583 p 3 
-21210055 p 4 + 27622695 p 5 - 279589 p 6 + 23497 p 7 ) 
2798298 + 103839446 p + 248730949 p 2 + 91081790 p 3 
-119321953 p 4 - 214483 p 5 + 2402 p 6 ) 
-987660 - 90012828 p - 421343630 p 2 - 300983008 p 3 
+354546288 p 4 + 19055623 p 5 - 2902607 p 6 + 156356 p 7 ) 
-164025 - 55522953 p - 52381 1627 p 2 - 718737692 p 3 
+737559653 p 4 + 137144016 p 5 - 24588235 p 6 + 1508228 p 7 ) 

1 1050744677 p 3 



+p 13 p 3 



+p 14 p 3 



+p 16 p 4 



+p 18 p 4 



+p 2 V 



+p 24 p 6 



331984 p 5 + 16780 p 6 , 
455796 - 3964687 p - 23412205 p 2 - 3663927 p 3 
+375501 p 5 + 2506 p 6 ^ 



1452867 p 4 



-23264478 - 483102731 p - 1256125498 p 2 + 
+534128347 p 4 - 99647028 p 5 + 3934654 p 6 + 455372 p 7 , 
-5993055 - 328891636 p - 1637004249 p 2 + 897069980 p 3 
+1365985483 p 4 - 226306853 p 5 - 13885678 p 6 + 5675792 p 7 ) 
-729000 - 161715291 p - 1608420692 p 2 + 82688597 p 3 + 2449767599 p 4 
-232309325 p 5 - 139531559 p 6 + 32366472 p 7 + 274104 p 8 ) 
-54779328 - 1192868328 p - 1036350248 p 2 + 3126514952 p 3 
+252124085 p 4 - 530576442 p 5 + 110376920 p 6 + 3523400 p 7 ) 
-11563290 - 661113157p - 1778966082 p 2 + 2728489572 p 3 
+1464568704 p 4 - 1261931324 p 5 + 243789711 p 6 + 21533484 p 7 ) 
-1166400 - 267060831 p - 1769360761 p 2 + 1310882671 p 3 
+3038379428 p 4 - 2084661610 p 5 + 338855326 p 6 + 82890944 p 7 ) 
-74824812 - 1221905422 p - 276377214 p 2 + 4066473373 p 3 
-2443698651 p 4 + 202709670 p 5 + 224730742 p 6 ) 
13126860 - 609332834 p - 1166508009 p 2 + 3871377708 p 3 



-1932786432 p 4 - 295096153 p 5 + 



t- 454945390 p 6 ) 
1103625 - 218082828 p - 1184923715 p 2 + 2650067734 p 3 



35 



-752045795/ - 1022894704 p 5 + 711726527 p 6 ) 
-P 25 P 7 (-53743284 - 754572613 p + 1226487190 p 2 + 416790514 p 3 

-1614449937 p 4 + 878716880 p 5 ) 
+p 26 p 7 (-8258625 - 336506619 p + 274671788 p 2 + 1010093115 p 3 

- 1758220080 p 4 +866894232 p 5 ) 
-p 27 p 7 (-605475 - 106456428 p - 91817943 p 2 + 978393477 p 3 

- 1446294368 p 4 +687560517 p 5 ) 
+p 28 p 8 (-23082921 - 120749456 p + 640575958 p 2 - 928058750 p 3 + 438708098 p 4 ) 
-P 29 P 8 (-3116070 - 61389864 p + 311935039 p 2 - 469461042 p 3 + 224128938 p 4 ) 
+p 30 p 8 (-200475 - 19494900 p + 116248678 p 2 - 186770845 p 3 + 90668832 p 4 ) 
-p 31 p 9 (-4042170 + 33260499 p - 57626524 p 2 + 28477028 p 3 ) 
+p 32 p 9 (-510300 + 7194033 p - 13394778 p 2 + 6717655 p 3 ) 
-3p 33 p 9 (-10125 + 376209 p - 740598 p 2 + 374614 p 3 ) 

+540p 34 (1 - p) p 10 (216 - 221 p) - 6075 p 35 (1 - pf p 10 . (A.2) 



A.3 V 3 (p,p,T) 



V 3 {p,P,T) 



2 + 2T + p (-6- 10T-4T 2 ) 



+pp 



5 + 12T -8T 3 +p 



-29 



+p 2 (78 + 138 T - 52 T 2 



, 2 2 
+P P 



+P 3 p 3 



52T + 24T 2 + 64T 3 + 16T 4 ) 
176 T 3 - 64 T 4 ) 
+p 3 (-36 - 60 T + 48 T 2 + 120 T 3 + 48 T 4 ) 
4T + 8T 2_ 8r 3 _ 16T 4 +p (_ 8 + 8T _24r 2 + 48T 4 ) 

+p 2 (56 - 22 T + 44 T 2 + 192 T 3 + 128 T 4 + 64 T 5 ) 
+p 3 (-262 - 336 T + 40 T 3 - 256 T 4 - 288 T 5 - 64 T 6 ) 
+p 4 (248 + 384 T - 120 T 2 - 384 T 3 + 128 T 4 + 384 T 5 + 128 T 6 

+p 5 (-66 - 110 T + 52 T 2 + 160 T 3 - 32 T 4 - 160 T 5 - 64 T 6 

- 10-40T-40T 2 + 8T 3 + 16T 4 

+p (58 + 164 T + 264 T 2 + 24 T 3 - 256 T 4 - 128 T 5 ) 
+p 2 (-206 - 378 T - 800 T 2 - 656 T 3 + 416 T 4 + 640 T 5 + 192 T 6 ) 
+p 3 (448 + 890 T + 1296 T 2 + 1080 T 3 - 256 T 4 - 864 T 5 - 384 T 6 ) 
+p 4 (-237 - 668 T - 1188 T 2 - 920 T 3 - 32 T 4 + 320 T 5 + 192 T 6 ) 
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+p 5 (-151 -70T + 760 T 2 + 904 T 3 + 256 T 4 + 32 T 5 ) 
+p 6 (162 + 230 T - 300 T 2 - 560 T 3 - 192 T 4 ) 

+p 7 (-36 - 60 T + 48 T 2 + 120 T 3 + 48 T 4 ) 

- 10 - 54T - 104T 2 - 56T 3 + 80T 4 + 64T 5 

+p (54 + 274 T + 492 T 2 + 400 T 3 - 144 T 4 - 416 T 5 - 192 T 6 ) 
+p 2 (-12 - 578 T - 1244 T 2 - 904 T 3 + 48 T 4 + 576 T 5 + 384 T 6 ) 
+p 3 (-206 + 456T + 2000 T 2 + 1688 T 3 + 304 T 4 - 160 T 5 - 192 T 6 ) 
+p 4 (156 - 356 T - 2060 T 2 - 2120 T 3 - 592 T 4 - 64 T 5 ) 
+p 5 (92 + 462 T + 1228 T 2 + 1272 T 3 + 400 T 4 ) 
-hp 6 (-88 - 270 T - 372 T 2 - 288 T 3 - 96 T 4 ) 
+p 7 (-22-2T + 28T 2 + 8T 3 ) 
+p 8 (30 + 46T + 16 T 2 ) 

+p 9 (-6- 10T-4T 2 ) 

- 3 - 20T - 56 T 2 - 64T 3 - 16T 4 + 64T 5 + 64T 6 

+p (11 + 160T + 432 T 2 + 376 T 3 + 80 T 4 - 96 T 5 - 128 T 6 ) 
+p 2 (-84 - 488 T - 1116 T 2 - 984 T 3 - 240 T 4 + 64 T 6 ) 
+p 3 (240 + 850 T + 1416 T 2 + 1224 T 3 + 336 T 4 + 32 T 5 ) 
+p 4 (-252 — 790 T - 960 T 2 - 712 T 3 - 208 T 4 ) 
+p 5 (74 + 318 T + 328 T 2 + 168 T 3 + 48 T 4 ) 
+/ (43 + 14T-28T 2 -8T 3 ) 
+p 7 (-33 - 48T- 16 T 2 ) 

+p 8 (6 + 10T + 4T 2 ) 1 (A.3) 



When T = — 1, we can further simplify the above expression 

V 3 (p, p, T = -1) = (1 - p) VP (l + p (1 - p) pf (l - p + p (1 - p) 3 ) 2 (A.4) 



+P 5 p 5 



There is no such factorization in the case T — 1: 



^ 3 (p,p J r = i) = 

4-20p + pp (9 + 23p-76p 2 + 120p 3 ) 

-2p 2 p 2 (6 - 12p - 231 p 2 + 583 p 3 - 384 p 4 + 110p 5 ) 

+p 3 p 3 (-66 + 126p - 792 p 2 + 2210p 3 - 2533 p 4 + 1731p 5 - 660p 6 + 120p 7 ) 
-2p 4 p 4 (40 - 234p + 865 p 2 - 1945 p 3 + 2518p 4 - 1727p 5 + 557p 6 - 6p 7 - 46 p 8 + 
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+p 5 p 5 (-31 + 835p - 2848 p 2 + 4098 p 3 - 2922/ + 936j> 5 + 21p 6 - 97 p 7 + 20p 8 ) (A.5) 
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